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We present a detailed study of the linear stability of plane Couette-Poiseuille flow in the 
presence of a cross-flow. The base flow is characterised by the cross flow Reynolds number, 
Ri n j and the dimensionless wall velocity, k. Squire's transformation may be applied 
to the linear stability equations and we therefore consider 2D (spanwise-independent) 
perturbations. Corresponding to each dimensionless wall velocity, k 6 [0, 1], two ranges 
of Ri n j exist where unconditional stability is observed. In the lower range of Ri n j, for 
modest k we have a stabilisation of long wavelengths leading to a cut-off Rinj- This 
lower cut-off results from skewing of the velocity profile away from a Poiseuille profile, 
shifting of the critical layers and the gradual decrease of energy production. Cross-flow 
stabilisation and Couette stabilisation appear to act via very similar mechanisms in this 
range, leading to the potential for robust compensatory design of flow stabilisation using 
either mechanism. 

As Ri n j is increased, we see first destabilisation and then stabilisation at very large 
Rinj. The instability is again a long wavelength mechanism. Analysis of the eigen- 
spectrum suggests the cause of instability is due to resonant interactions of Tollmien- 
Schlichting waves. A linear energy analysis reveals that in this range the Reynolds stress 
becomes amplified, the critical layer is irrelevant and viscous dissipation is completely 
dominated by the energy production/negation, which approximately balances at critical- 
ity. The stabilisation at very large R{ n j appears to be due to decay in energy production, 
which diminishes like R^j- ^ ur study is limited to two dimensional, spanwise indepen- 
dent perturbations. 



1. Introduction 

From the perspective of applications in technology, Poiseuille flow of viscous fluid along 
a duct is undoubtedly one of the most important flows studied as it underpins the field 
of hydraulics. Instability and subsequent transition from laminar flow marks a paradigm 
shift in the dominant transport mechanisms of mass, momentum and heat, and it is for 
this reason that the subject remains of enduring interest, even after more than 100 years 
of study. In this paper we focus on two methods for affecting the linear stability of plane 
Poiseuille (PP) flow. The first method consists of introducing a Couette component to 
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the flow, by translation of one of the walls. The second method consists of introducing a 
cross-flow, e.g. via injection through a porous wall. While both effects have been studied 
individually to some extent, there are fewer studies of the two effects combined, which is 
the main focus here. 

We first summarise the effects of a Couette component on a plane Poiseuille flow. 
The main curiosity here stems from the observation that PP flow is linearly unstable 
when the critical Reynolds number exceeds R c ~ 5772, (Reynolds number based on the 
center-line axial velocity and the half- width of the channel; see |Orszag (1971)) , whereas 
the plane-Couette (PC) flow is absolutely stable with respect to infinitesimal amplitude 
disturbances, R c — oo; see Romanov (1973)[ Superimposing PP and PC flows, we may 



ask if a small Couette component can affect the stability of the PP flow. Stability of 
plane Couette-Poiseuille (PCP) flow was first studied by Potter (1966)| and later by 
Hains (1967)] |Reynolds k Potter (1967)1 and |Cowley & Smith (1985)] The results are 



typically understood with respect to a Reynolds number that is based on the maximal 
velocity of the Poiseuille component, say R p , and the ratio of wall velocity to maximal 
velocity of the Poiseuille component, denoted k. For small Couette components, k, it is 
possible to observe some destabilisation of the flow, (depending on the wavenumber), but 
as soon as k > 0.3 a strong stabilisation of the flow sets in. As the velocity ratio k exceeds 
0.7, the neutral stability curve completely vanishes and the flow becomes unconditionally 
linearly stable, i.e. R c — > oo. The term "cut-off" velocity has been used to describe this 



stabilisation; see Reynolds & Potter (1967) 



Although plane Couette flows are widely studied it is worth noting that they are ac- 
tually difficult to produce, i.e. outside of the computational and theoretical domain. In 
many duct flows axial translation of a wall is either not possible or is limited in terms of 
speed. High R frequently means high velocities, lowering the range of achievable k as the 
flow velocity increases. Therefore, the range of practical flows for which a sufficiently sta- 
bilizing Couette component can be introduced is limited and we know of no technological 
applications where this is used for stabilisation. 

Annular Couette-Poiseuille (ACP) flows are more practically relevant and have also 
been studied extensively (Mott & Joseph (1968) , Sade ghi fc Higgins (1991)| ). For example, 
ACP flows occur when removing/inserting drillpipe or casing from vertical wellbores 
during an operation called "tripping" . Sadeghi and Higgins studied the flow between two 
concentric cylinders, the outer being stationary while the inner is moved with a constant 
(dimensionless) velocity k in the streamwise direction. They showed that varying the 
radius ratio (?]) between the outer and inner cylinders can have a dramatic effect on the 
stability characteristics. The limit 77 — ?> 1 approximates PCP flow and is unconditionally 
linearly stable for k ^ 0.7, thereby confirming |Potter (1966)[ By increasing 77, the cut-off 
condition is attained for lower values of k and the cut-off relation between k-r\ is almost 
linear. Similar to Mo tt fc Joseph (1968) [ they argued that increasing 77 increases the 
asymmetry of the base flow profile which in turn increases the stability with respect to 
axisymmetric disturbances. Their findings are very relevant to our work, since we later 
show that the stability achieved by increasing 77 in ACP flows and that achieved by 
applying a small cross-flow in PCP flows are essentially similar. 

Shear flows with cross-flow occur in a range of natural settings as well as in vari- 
ous technological applications. As examples, we cite studies in sediment-water interfaces 
over permeable seabeds ( |Goharzadeh et al. (2005) ); fluid transport and consequent mass 
transfer at the walls of blood vessels, the lungs and kidneys (Majdalani et al. (2002) ); and 
flow through the fractures of geological formations ( |Berkowitz (2002) \ . In some techno- 
logical applications cross-flow is an inherent part of the process, e.g. dewatering of pulp 
suspensions in paper making, whereas in others it is introduced to affect the stability. 
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An example of the latter is the use of wall suction to delay the transition to turbulence 
over the surface of an aircraft wing (Joslin (1998)). 

The stability of PP flow with cross- flow was first analysed by |Hains ( 1971) | and Shep- 



|pard (1972)| both of whom have shown that a modest amount of cross-flow produces 
significant increase of the critical Reynolds number. These results are however slightly 
problematic to interpret in absolute terms, since at a fixed pressure gradient along the 
channel, increasing the cross-flow decreases the velocity along the channel, (hence effec- 
tively the Reynolds number) . This difficulty was noted by |Fransson fe Alfredsson (2003) [ 
who used the maximal channel velocity as their velocity scale (instead of that based on 
the PP flow without cross-flow), and thus separated the effects of base velocity magnitude 
from those of the base velocity distribution. Using this velocity scale in their Reynolds 
number R, they showed regimes of both stabilisation and destabilisation as the cross- 
flow Reynolds number was increased. For example, for R = 6000 and wavenumber a = 1, 



Fransson fc Alfredsson (2003)| have shown that the cross-flow was stabilising up to a cross- 
flow Reynolds number R in j « 3.4, and then starts destabilising before re-stabilising again 
at Ri n j ~ 635. The initial regime of stabilisation is the one corresponding to the earlier 
results. 

Although cross-flow affects the base velocity profile, the main change to the linear 
stability problem is to add an inertial cross-flow term to the Orr-Sommerfeld operator. 
One reason why addition of terms like the cross-flow term can destabilise an otherwise 
stable shear flow is suggested by the two-dimensional instability of Blasius boundary 
layer, as studied by |Baines et al. (19 96) In such flows, the resulting growing disturbance 
is known as a Tollmien-Schlichting (T-S) wave. They showed that the interaction between 
two idealized modes, viz. an "inviscid" neutral mode at zero viscosity and a decaying 
viscous mode (or modes) existing at uniform shear undergo resonant interactions. The 
latter is forced by the former through the no-slip wall boundary conditions. 

In the present study we focus on the combination of cross-flow and Couette component. 
Our motivation stems from a desire to understand how the two mechanisms interact, since 
in terms of technological application different mechanical configurations may be more or 
less amenable to cross-flow and/or wall motion. This means that there is value in knowing 
when one effect may compensate for the other in stabilizing (or destabilizing) a given 
flow. To the best of our knowledge, stability of PCP flow with cross-flow has only been 
studied in any generality by |Hains (1971)| In considering the base flow for PCP flow 
with cross-flow, (which is parameterized by Ri n j and k), the relation kRi n j — 4, defines 
an interesting paradigm in which the base velocity in the axial direction is linear. These 
Couette-like flows have been studied by Nicoud & Angilella (1997) for increasing Ri n j. 
They found a critical value of Ri n j ~ 24, below which which no instability occurs, (we 
have translated their critical value of 48 into the Ri n j that we use). Therefore, we observe 
that understanding of cross-flow PCP flows is far from complete. We aim to contribute 
to this understanding. 

The 3D linear stability of PCP flow with cross-flow is amenable to Squires transfor- 
mation, so that the linear instability occurs first for 2D (spanwise-independent) pertur- 
bations. It is these perturbations that we study here. Our aim is to demarcate clearly 
in the (Rinj, /c)-plane, regions of unconditional stability, i.e. where there is a cut-off wall 
velocity or injection velocity. We also wish to understand the underlying linear stability 
mechanisms as Ri n j and k are varied. 

Although the study of 2D perturbations is justified from the pure perspective of linear 
stability, it must be acknowledged that 3D and nonlinear effects are likely to be relevant 
in instabilities that are observed to grow, i.e. the actual transition. The past two decades 
have seen extensive study of transient growth mechanisms, due to non-normality of the 
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operator associated with linearized Navier-Stokes equations. Algebraic growth of 0(Re 2 ) 
may occur for linearly stable disturbances that decay only slowly over a timescale of 
O(Re). It has been proposed that this transient algebraic growth is responsible for sub- 
critical transition in wall-bounded shear flows. For an overview of these developments 
we refe r to |Reddy et al. (1993)|[Schmid fc Henningson (2001)[ |Chapman (2002)[ |Schmid| 



(2007)[ 



At the same time as transient growth mechanisms have undergone extensive research, 
self sustaining nonlinear mechanisms were proposed by Waleffe and others, e.g. Hamilton 
et al. (1995) , Waleffe (1997) | In this scenario energy from the mean flow can be fed 



back into streamwise vortices, thus resisting viscous decay. Self-sustained exact unstable 
solutions to the Navier-Stokes equations were found by Fai sst fc Eckhardt (2003)] and 



by Wedin k, Ker swell (2004)[ Much current effort is focused at understanding the link 



between these self-sustained unstable solutions and observed transitional phenomena, 
such as intermittency, streaks, puffs and slugs; see e.g. Hof et al. (2004)||Hof et al. (2005) [ 
Eckhardt et al. (2007)[ |kerswell fc Tutty (2007)| 



Our study does not deal with any of the complexities of transition mentioned above, 
and as such the relevance may be questioned. This is a fair criticism, but on the other 
hand we note that for other classical shear flows that are linearly stable at all Re, careful 
control of apparatus imperfections and the level of flow perturbations can significantly 
retard the point at which transition is observed. For example, in Hagen-Poiseuille flow of 
Newtonian fluids one typically observes transition to turbulence starting for Re > 2000. 
However, an experimental flow loop in Manchester UK produces stable laminar flows for 
Re » 24, 000, |Hof et al. (2004)[ |Peixinho fe Mullin (2006)] and stable flows have even 



been reported up to Re ~ 100, 000, Pfenniger (1961) This all suggests that significantly 



enhanced stability may be achieved experimentally, where predicted by the linear theory. 

The question of how to achieve a PCP flow with cross-flow in practice is also rele- 
vant. Evidently all Poiseuille flows occur in finite geometries with entrance effects, side 
walls and imperfections in the planar walls, so that the notion of a truly planar infi- 
nite flow is anyway flawed. Uniform base flows studied in hydrodynamic stability are 
invariably an approximation of experimental reality. Even in the absence of wall motion 
planar Poiseuille flow is difficult experimentally, due to spanwise perturbations and in- 
flow non-uniformities. This said, a geometry with a uniformly translating channel walls 
is particularly difficult to achieve and as mentioned before, k ~ 0.7 is difficult for high 
R flows where R is increased via flow rate. Imposing a uniform cross-flow along with a 
streamwise pressure variation is more practically achievable in practice. As an example, 



see|Vadi fc Rizvi (2001) A uniform trans-membrane pressure cross-flow micro-filtration 



system is able to maintain uniform trans-membrane pressure with high cross-flow velocity 
(Yinj) an d improves the utilization of available filtration area. In the patent of 



Sandblom 



(2001)| the concept of operating a membrane filtration unit using UTMP has been pro 



posed, such that pressure drop along the channel can be adjusted independent of the 
Vi n j. A different generic concept for achieving a uniform cross-flow over a finite length of 
a porous-walled channel, consists of injecting the fluid along a secondary channel behind 
one of the porous wall that has a linearly converging geometry. Although there are clear 
technical challenges, it is worth remarking that the injection velocities needed for sta- 
bility are very modest by comparison with the wall velocities, and therefore as a target 
appear achievable. 

An outline of our paper is as follows. We commence below in ^2] by introducing the 
base flow and linear stability problem. We describe our numerical method and present 
benchmark results that illustrate typical effects of varying Ri n j and k. These results serve 
to motivate the presentation of results, which follows in the following 3 sections. Section[3] 
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considers low Ri n j and significant k, where we see that long wavelengths dominate. In £j4] 
we characterise the flows for intermediate Ri n j and small k, where short wavelengths are 
the least stable. Finally, we consider large Ri n j in $5 where we find both destabilisation 
and eventual stabilisation. The paper concludes in £6] with a summary of the principal 
results. 



2. Stability of plane Couette-Poiseuille flow with cross-flow 

The base flow considered in this paper is a plane Couette-Poiseuille flow (PCP) with 
imposed uniform cross-flow. This flow is two-dimensional, viscous, incompressible and 
fully developed in the streamwise direction, x,(all dimensional variables are denoted with 
a "hat", i.e. : ). The imposed base velocity in the y-direction, v, is constant and equal to 
the injection/suction velocity V{ n j. Since v is constant, the x-component of velocity, u 
depends only on y. The flow domain is bounded by walls at y — ±h, and is driven in 
the indirection by a constant pressure gradient and by translation of the upper wall, at 
speed U c . The x-component of velocity, u(y), is found from the ^-momentum equation, 
which simplifies to: 

du 1 dp d 2 u 



oy p ox oy z 



(2.1) 



where p is the density, v — p/p is the kinematic viscosity, and p the dynamic viscosity. 
The boundary conditions at y — ±/i are: 



l(-h) = 0, u(h) = U c 



(2.2) 



To scale the problem we scale all lengths with h, hence (x,y) — (x/h,y/h). For the 
velocity scale two choices are common. First, the imposed pressure gradient defines a 
"Poiseuille" velocity scale: 

h 2 dp 
2p dx ' 

which is equivalent to the maximum velocity of the plane Poiseuille flow, driven by the 
pressure gradient alone. Secon d, we may tak e the maximum velocity, which we need to 
compute. Up is the choice of Potter (1966) and thus allows one to compare directly 



(2.3) 



with the studies of PCP flows. In the absence of a cross-flow, the maximal velocity is not 
actually very sensitive to the wall velocity U c , at least for U c < U p , which covers the range 
over which the flow stabilises. However, in the case of a strong cross-flow the x-velocity 
is reduced significantly below U p , which therefore loses its meaning. Consequently we 
adopt the second choice and scale with the maximal velocity, U max . This choice retains 
physical meaning in the base velocity, but does introduce algebraic complexity. 



The solution is found from (2.1)-(2.2) after detailed but straightforward algebra: 



AcoshR inj -kR inj e- R ^ + [kR inj -4\e R ^y/ h 2y 



where k and Ri n j are defined by: 



2Rinj sinh Rj jTl j 



U p 

V- h 



(2.4) 

(2.5) 
(2.6) 
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These two dimensionless parameters uniquely define the dimensionless base flow. The 
parameter k is the velocity ratio of Couette to Poiseuille velocities, which is useful as it 
allows direct comparison with earlier results on stabilisation of PCP flows without cross- 
flow. The parameter Ri n j is simply a Reynolds number based on the injection velocity. 
Primarily here we consider the ranges: k E [0, 1] and Ri n j ^ 0. 

For relatively weak crossflow velocities, the velocity component u(y) has a single max- 
imum at a value of y — y m ax defined by: 



sinhi? 



inj 



Ri 



kR, 



(2.7) 



The maximal velocity U max , is then evaluated from (2.4). Since, smhRi n j ^ Rinj, we 
can see that y max > for k ^ and Ri n j > 0. Both the injection cross-flow and Couette 
component act to skew the velocity profile towards the upper wall. For stronger cross- 
flow velocities, (or sufficiently large k), the maximal velocity occurs at the upper wall, 
i.e. Umax — U c . 

The division between weak and strong cross-flows, taking into account also the Couette 
component, is defined by the line 



kRinj 4 



_ sinh R inj 

Rinj&^ tn ^ 



(2.8) 



The dimensionless base velocity is given by: 



u(y) = < 



4 cosh R 



tv j 



kR inj e- Rin J - [4 - kR i nj]e RinjV + 4y sinhii 



4 cosh Ri nj - kR inj e R ^ - [4 - kRi n j]e RinjVma ° + 4y max sinh^ nj 

sinh i?,-„ 



4 cosh Ri n j — kRi n je~ 



[4 - kR in j]e 



kRinj ^ 4 



RinjV 



■Rime 

2y 



^Rinj sinh Rinj 



R 



(2.9) 



kRinj ^ 4 



1 



sinh Ri 



It can be verified that in the limit R 
velocity profile is retrieved: 



I 71 J 



R inj e n » 

0, with k fixed, the classical form of PCP base 



l-y 2 + ±{l + y) + R mj [\{y 



hi 

16 



Rinj [ q 



61- 



(2.10) 



as R^j -> 0, with k < 4[1 - sinh. R inj / {Ri n jQ Rtn3 )]/ Rinj ~ 4[1 + R in j/3]- 

Examples of the base velocity profile are given in Fig. [l] for for k — 0.5 and different 
values of Rinj- Observe that for Ri n j — 8, when kRinj — 4, the velocity profile is linear. 
This flow has been termed a "generalised Couette" flow by |Nicoud h Angilella (1997)| 

We shall denote differentiation with respect to y by the operator D. The first and 
second derivatives of the base flow, Du and D 2 u respectively, influence the stability of 
the flow. We find that D 2 u has sign determined by (kRi n j —4), and increases in absolute 
value exponentially towards the upper wall. For kRi n j < 4, the velocity is concave, and 
is convex otherwise. Since D 2 u does not change sign, the maximal absolute value of the 
first derivative is found either at the upper or lower wall, y — ±1. The maximal velocity 
gradients are found at the lower wall for small Rinj and also for a range of Rinj close to 
kRi n j — 4, but otherwise are found at y = 1; see Fig. [2|a). At large Ri n j the maximal 
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Figure 1: Mean Velocity Distribution for k — 0.5 and Ri n j— 0, 1, 4, 8, 15 and 30 (Ri n j — 
30 marked with a □) 




Figure 2: (a) Maximal velocity gradient, \Du\ max , plotted against Ri n j for k — 
0.35, 0.5, 0.65, (k — 0.65 marked with a □). The thick line indicates where the max- 
imum is attained at y = — 1; otherwise at y = 1. (b) Variation of _D 2 w with y for k — 0.5 
for : # inj = 0, (□); R inj = 4, (o); R inj = 8, (x); i^ nj = 12, (o). 



velocity increases almost linearly: 

~ Rinj — 



[kR inj - 4]e J 



2 sinh i?,- 



fc kR 



0(e 



')■ 



7?; 



(2.11; 



Figure [2^b) shows examples of the profiles of D 2 u. We observe that -D 2 u ~ over a 
large range of y, close to the lower wall, whenever a significant amount of cross-flow is 
present, i.e. R inj > 1. 
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2.1. Dimensionless groups 

The base base flow is fully defined by the parameters k and Ri n j, as discussed above. In 
addition, the transient flow and associated stability problem will depend on the stream- 
wise Reynolds number, R, which we define in terms of U max , i-e. 

U h 

*J m,a:r, ' » 



R 



(2.12) 



To aid the reader in interpreting our results in terms of those previously published, it 
is helpful to consider also a Reynolds number based on the Poiseuille velocity, U p , say 
Rp — Uph/u. Straightforwardly, we find R — R p F(k, R{ n j): 



4 cosh R 



> n j 



kR 



inj 5 



+ [kR inj - 4]e 



RinjVr 



4y r 



sinh R 



inj 



F{k,R inj ) = < 



^Rinj Sinh Ri n j 
kRinj ^ 4 

kRinj ^ 4 



sinh Ri 



sinh Rj 



, (2,13) 

Note that F(k, Rinj) — U max /U p , which is fixed by the parameters k and Ri n j- Thus, for 
fixed k and Ri n j an increase in R is interpreted as an increase in R p , and vice versa. It is 
also useful to know the ratio of upper wall velocity to the maximal velocity, i.e. U c /U max , 
which we shall denote k, given simply by the ratio k/F(k, Ri n j)'- 



2k Ri n j sinh Ri n j 



4 cosh Ri 



kRi 



k(k, Rinj) — < 



[kRir, 
kRinj ^ 4 

kRinj ^ 4 



^ sinh R 
1 



sinh Ri 



The ratio R/R p and the upper wall speed k are illustrated in Fig. [i]for 



sinh Ri 



(2.14) 
convenience. 
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2.2. The stability problem 

The base flow is two-dimensional, but since v — Rinj is constant, the 3D linear stability 
equations are only modified by the addition of a constant convective term: 

R *Lu' 

where u' = (u',v',w f ) denotes the linear perturbation. The classical Squire transforma- 
tion can therefore be applied to the temporal problem, showing that for any unstable 3D 
linear disturbance there exists an unstable 2D linear disturbance at lower R; see Squire 
(1933)[ 

It suffices to consider only 2D disturbances and we adopt the usual normal mode 
approach to linear spatially periodic perturbations, introducing a stream function which 
we represent in modal form as: 

j,(x,y,t)=(i>(y)e^- ct )\ (2.15) 

with u' — D(j)(y)e^ ia ^ x ~ ctS)li , v' — —ia(f)(y)e^ a ^ x ~ ct ^ . Thus, a is real, denoting the wavenum- 
ber, c denotes the complex wave speed, (c = \/— T), and (j)(y) denotes 

the amplitude of the stream function perturbation. The modified Orr-Sommerfeld (O-S) 
equation for the flow is: 

u)(a 2 - D 2 ) - D 2 u](j) - R in jD(a 2 - D 2 )cj) = (a 2 - D 2 ) 2 (j), 



iaR[(c 

and the boundary conditions are 



;±i) = D0(±i) = o. 



(2.16) 



(2.17) 



The inclusion of the injection cross-flow results in an additional 3rd order derivatives 
in the inertial terms, i.e. Ri n jD(a 2 — D 2 )(j). Note that Ri n j also influences stability via 



the base velocity profile u(y). Equations (2.16 )-(2.17) constitute the eigenvalue problem. 
The eigenvalue c is parameterised by the 4 dimensionless groups (a, R, Ri n j, k) and the 
condition of marginal stability is: 



Ci(a,R,R inj ,k) = 



(2.18) 



We attempt to characterise the stability of (2.16 )-(2.17) for positive (a, R, Ri n j 7 k). We 



may note that the limit R — > oo for finite Ri n j, reduces (2.16) to the Rayleigh equation. 
Since D 2 u is of one sign only, there are no inflection points and hence no purely inviscid 



instability. This suggests that the instabilities of ( 2.16 )-( 2.17) will be viscous in nature 



Addition of the constant cross-flow terms does not fundamentally alter the O-S prob- 



lem, and we expect a discrete spectrum. To find the spectrum of (2.16)-(2.17) we use a 
spectral approach, representing by a truncated sum of Chebyshev polynomials: 



J2 a nT n (y) for y E [-1,1] 



(2.19) 



n=0 



where TV is the order of the truncated polynomial, a n is the coefficient of the n-th Cheby- 
shev polynomial, T n (y). This method is described for example in Schmid & Henningson 



(2001) and is widely used. The discretised problem is coded and solved in Matlab. The 
accuracy of the code has been checked against the results of Mack (1976) for the Bla- 
sius boundary layer, with various results for PP flow in |Schmid fc Henningson (2001)[ 
with the PCP flow results of |Potter (1966) and finally against results for PP flow with 
cross-flow; see |Sheppard (1972)[ |Fransson fc Alfredsson (2003) | Th e results are accurate 
up to three, four and five significant places when validated against |Potter (1966)| |Mack 
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|(1976)| and |Fransson fc Alfredsson (2003) | respectively. All the numerical results given 
below have been computed with N — 120. On using 200 collocation points, the growth 
rates changed only in the fourth significant place in the worst case. 



2.3. Characteristic effects of varying k and R 



> n j 



Before starting a systematic analysis of (2.16 )-(2.17), we briefly show some example re- 
sults that illustrate the characteristic effects of varying the Couette component, k, and 
the cross-flow component, Ri n j- These examples also serve to establish the framework 



of analysis used later in the paper. With reference to PP flow, Potter (1966) first ob- 
served that the stability is increased by adding a Couette component while |Fransson 



Alfredsson (2003) | showed that cross-flow can stabilise or destabilise PP flow. 



2.3.1. Eigenspectra 

Setting (a,R) — (1,6000), we investigate variations in the eigenspectrum of (2.16)- 



(2.17). According to a classification proposed by Mack (1976) the spectrum of PP flow 
spectra may be divided into 3 distinct families: A, P and S. Family A exhibits low 
phase velocity and corresponds to the modes concentrated near the fixed walls. Family 
P represents phase velocities, c r , close to the maximum velocity in the channel. Family 
S corresponds to the mean modes and has phase velocity c r close to the mean velocity. 
In Figsga) & (b), we track the eigenmodes as k and Ri n j, respectively, are varied from 
zero. The initial condition (denoted by □) represents the PP flow. 

Referring to Fig. ga), (where Ri n j — 0), addition of the Couette component increases 
the mean velocity: the S modes shift from c r — 0.6667 at k — (PP flow) to c r — 0.7513 at 
k = 1. The family of P modes is also shifted to the right. The A modes are associated with 
both walls and as k increases we see a splitting of the family, with the upper wall modes 
moving to the right as k is increased. The least stable mode is a wall mode associated with 
the lower wall, which we observe stabilises monotonically as k is increased. Figure gb) 
shows the effects of increasing Ri n j, (holding k — 0). The least stable A mode of PP flow 
initially stabilises and then destabilises with increasing Ri n j. This behavior has also been 



observed by Fransson & Alfredsson (2003) For large Ri n j the A, P, and S families have 



disappeared, instead leaving two distinct families of modes. It appears that each of the A, 
P, and S families splits, with some modes entering each of the two families (this alternate 



splitting is most evident for the S modes). As observed by|Nicoud fc Angilella (1997) 



the phase speed no longer lies in the range of the axial velocity. This does not violate 
the conditions on c r , given by Joseph (1968) and Joseph (1969), since these conditions 
are derived for parallel flows only. 



2.3.2. Increasing R 



i ii j 

Next, we illustrate the qualitative effects of increasing R in j at fixed (a,R,k), in 
Fig. ga). We again fix a — 1 and R — 6000, and show the variation of the least stable 
eigenvalue, for k — 0, 0.5, 1. Our results for k — (PP flow) may be compared directly 
with those of Fransson & Alfredsson (2003) We observe that as Ri n j increases we have 
an initial range of stabilisation (cj jCr jt decreasing), followed by a range of destabilisa- 
tion (ci jCr it increasing), and finally again stabilisation at large Rinj, (c^crit decreasing). 
Qualitatively, we have observed these same three ranges of decreasing/increasing c^nt, 
as Ri n j increases, for all numerical results that we have computed, and this provides a 
convenient framework within which to describe our results. 

For fixed (a, R, k), the case Ri n j = may either be stable or unstable, in which cases 
there are respectively two or three marginal stability values of Rinj- We denote these 
marginal values of Ri n j by: Ri n j,i, Rinj,2, Rinj, 3, noting that in the case that Ri n j — is 
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Figure 4: Eigenspectrum of (a,R) — (1,6000) by varying k and Ri n j- 40 least stable 
modes are considered, (a) Effect of increasing k from to 1 in steps of 0.01, keeping 
Rinj — 0. (b) Effect of increasing Ri n j from to 100 in steps of 0.05, keeping k = (PP 
flow). The symbols in (a) and (b) are similar and are denoted as follows: k — or R{ n j — 
by (□), k — 1 or Ri n j — 100 by (o) and intermediate k or R{ n j by (.). Note that the 
PP flow spectrum is represented by the □ in both figures, and shows the vertical family 
of S-modes, the branch of A-modes (diagonally upwards from centre to left) and branch 
of P-modes (diagonally upwards from centre to right). 
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Figure 5: (a) Effect of increasing Rinj on the stability of PCP flow, for (a, R) — (1, 6000) 
and different values of k — (□), 0.5 (o), 1 (x). (b) Maximal growth rate for increasing 
Rinj at different R, (k — 0.5 and the step in values of R between curves is 10 4 ). 



stable Rinj,! is absent. More clearly, Ri n j,2 will always represent a transition from stable 
to unstable, while Ri n j,i & Rinj,3 denote transitions from unstable to stable. The PCP 
flows for k — 0.5 and 1 are stable for (a,R) — (1,6000) in the absence of cross-flow, 
Rinj — 0. For a larger R, k — 0.5 is unstable at Ri n j — 0, but k = 1 remains stable for 
all (a, R). 

Figure [5|b) shows the maximal growth rate 7, for increasing Ri n j at different R, with 
k — 0.5. The maximal growth rate is computed over wavenumbers a £ [0, 1]: 

7 = max {aci}, (2.20) 

«G[0,1] 

which often captures the largest growth rates over all a. We observe that the first marginal 
value Rinj, 1 increases with R, but appears to converge towards a finite value as R — > 00. 
The second marginal value of Ri n j,2 appears independent of R, (at least numerically). 
For k — 0.5 we have Ri n j,2 ~ 24.7. Nicoud fc Angilella (1997)] have observed a similar 
behaviour in studying the generalised Couette flow, (for which the constraint, kRi n j — 4, 
is always satisfied). They have found Ri n j,2 ~ 24 (note that Nicoud and Angilella use 
the full channel width as their length-scale, and therefore report Ri n j,2 ~ 48, in their 
variables). In contrast, the third marginal value, Rinj, 3, is strongly dependent on R. For 
example, for k — 0.5, the values corresponding to R = 10000 and 100000 are Ri n j,3 ~ 83 
and Rinj,3 ~ 287 respectively. 

2.3.3. Increasing k 

Figure [6] explores the effects of increasing the Couette component k, on 7 and on the 
marginal values of Rinj- Figure[6|a) indicates that the sensitivity of Ri n j,2 to k is also not 
extreme: we have found that this transition occurs within the range ~ 22—25 for k £ [0, 1]. 
For each value of k examined, we also observe numerically a similar independence of 
Rinj,2 to R as seen earlier in Fig. [5^b) for k — 0.5. The 3rd marginal value, Ri n j,3, 
is strongly dependent on k. For example, at R — 40000, Ri n j,3(k — 1) ~ 120 and 
Rinj,3(k — 0.8) ~ 135. In general, increasing k shifts Ri n j,3 to the left, thereby decreasing 
the span of the unstable region. Increasing k also decreases the maximum value of 7. 
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Figure[6jb) looks at the first transition, Ri n j,i at R = 40000. Potter (1966) was the first 
to observe that for PCP flows (i.e. Ri n j — 0), a gradual increase in the wall velocity results 
in crossing a "cut-off" value of k, say k\, such that for k > k\ the flow is unconditionally 
linearly stable. It has already been pointed out from the results of Fig. [5|b)that Ri n j,i 
is finite as R — > oo. In addition, the results in Fig. [6|b) indicate that Ri n j,i decreases 
with A; at a finite R. Hence, it can be inferred that as R — > oo, the cut-off wall velocity, 
ki — ki(Ri n j) must decrease with Ri n j. 



3. PCP flows and the effects of small Ri n j 

Having developed a broad picture of the different transitions occurring in the flow, we 
now focus in depth at each range of Rinj, to understand the stability mechanisms in play. 
We start with the range of small Ri n j- 

PCP flows without cross-flow are stable to inviscid modes, but viscosity admits addi- 
tional modes, i.e. the Tollmien-Schlichting (TS) waves, which may destabilise, according 
to the value of k. When aR 3> 1 with c ~ 0(1), viscous effects occur in thin oscillatory 
layers: (i) adjacent to the walls, (of thickness ~ (ai?) -1 / 2 ), and (ii) close to the critical 
point(s), y C7 where u(y c ) — c r ^ cr i t are found, (of thickness ~ (ai?) -1 / 3 ). It is in the crit- 
ical layers that we see peaks in the distribution of energy production, implying transfer 
from the base flow. |Potter (1966)] put forward the argument that for a dimensionless wall 
velocity that exceeds c r;Cr ^, the critical layer near the moving wall will vanish and there 
remains only one critical layer, near the fixed lower wall. The thickness of this second 
layer increases with wall velocity, thereby favouring stabilisation. 

This mechanism appears to correctly describe the long wavelength perturbations, (at 



Rinj = 0), which are found to be the least stable for k ~ 0(1). Indeed Cowley & 
Smith (1985) developed a long wavelength analysis (a ~ R) 7 in order to predict the 
cut-off value k\(Ri n j — 0) ~ 0.7. For values k ~ 0(1), PCP flows have only a single 
neutral stability curve (NSC). However, Cowley Smith ( 1985) | noted that for smaller 
k, multiple neutral stability curves could exist, and at shorter wavelengths. For example, 
when ^ k ^ R~ 2 / 7 there is one NSC, when R~ 2 / 7 ^ k ^ R~ 2 / 13 there are thr ee NSC's, 
and when R~ 2 / 13 ^ k < 1 there are two NSC's; see|Cowley fc Smith (1985)[ Thus to 
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understand the effect of cross-flow in PCP flows, the different regimes of k need to be 
considered separately. 

For R in j « 0, we expect the stability behaviour to be close to that of the PCP flow 
without cross-flow. Intuitively we expect cross-flow to stabilise, and so study the range 
c r ,crit < k ^ k\(Ri n j — 0). We examine the NSC's obtained from the O-S equation 
corresponding to k = 0.5, under different values of Ri n j] see Fig. [7^a). As expected, 
increasing R in j results in a progressively larger critical R — R cr it- We also observe 
that both the upper and the lower branches are oriented at an angle of 45 degrees, 
(i.e. a ~ -R -1 ), at high values of R. On fixing Ri n j and increasing k we have found that 
for successively large k the upper and lower branches move together as R cr it increases, 
eventually coalescing at k = k\(Ri n j). This mechanism is identical with that observed by 
Cowley fc Smith (1985)| suggesting the applicability of a long wavelength approximation 



in order to predict ki(Ri n j). Figure [7|b) plots the values of c r at criticality, as Ri n j is 
varied, also for k — 0.5. The critical values are tabulated in Table [I] The dependence is 
initially linear. We observe that k > c rjCr a over the computed range. 

3.1. Long wavelength approximation 



We follow the long wavelength distinguished limit approach of Cowley & Smith (1985) 
taking a — > and R — > oo, with A = (ai?) -1 fixed. The product aR is fixed along the 
upper and lower branches of the NSC. Thus, as the two branches of the NSC coalesce 



in the (k 7 A)-plane we observe k — > ki(R in j). In the long wavelength limit, equation 2.16 
becomes: 

iX [D 4 - R inj D 3 ] cj) + {u-c) D 2 (j) - {D 2 u)(j) = 0, (3.1) 
with boundary conditions ( |2.17 ). 



Figure[8|a) shows the NSC obtained from (3.1), plotted in the (k, A)-plane for various 
Rinj- The cut-off value k\(Ri n j) is the maximal value of k on each of these curves. 
These values are listed in Table [2] We also list the dimensionless wall speeds at cut-off, 
i.e. k(k\, Ri n j). We observe that the cut-off wall speed decreases with Rinj- This is in 



agreement with the concluding remarks of ^2.3.3 
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Table 1: Critical values for k — 0.5 and increasing Ri 
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Figure 8: (a) Long wave NSC's showing the dependence of A on k for R in j —0 (dashed 
line), 0.3 (dash-dot), 0.5 (solid), 0.7 (dash-dot-dot) and 1 (long dash); (b) Ci ;Cr i t versus k 
for A = 2.5 x 10~ 5 , and R in j =0 (dashed line), 0.3 (dash-dot), 1 (solid), 1.2 (dash-dot-dot) 
and 1.3 (long dash). 



Figure |8[b) shows Cj for the least stable eigenvalue of the long wavelength problem, 
for fixed A = 2.5 x 10~ 5 and different values of Ri n j, as k is varied. When R in j ^ 1.3, we 
find that Ci^ cr n ^0, V k G (c rjCr ^, fei(0)], implying that there are no neutral or unstable 
long wavelength perturbations in this range of k, (i.e. at least until we approach the 
second transition at Rinj, 2)- Thus, in this initial range of say Ri n j ^5 1-3, provided that 
k > c rjCr it, we can talk equally of a cut-off value for k or for Ri n j- 
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Table 2: Cut-off values, k\, and wavespeed c rjCr i t , for increasing Ri 



3.2. Effects of asymmetry of the velocity profile 

We observe that Ri n j enters the stability problem in two distinct ways. The first one 
represents the direct contrib ution of the additional third order inertial term, Ri n jD(a 2 — 
D 2 )(j), in the O-S equation (2.16). For the second one, Ri n j influences the base velocity 
pr ofile. T o exp lore which of these effects is dominant, we show in Fig. [9] the spectra 
of (2.16 )-(2.17) obtained with and without the term, Ri n jD(a 2 — D 2 )(f), included in the 
computation. The critical parameters corresponding to Ri n j — 0.5, in Table[TJ are chosen 
and fixed for this comparison. Figure[9|a) shows the two spectra at Ri n j — 0.5, which are 
near identical, completely overlapping on the figure. This suggests that at smaller Ri n j, 
the effects of cross-flow manifest completely via the base flow velocity profile. Figure 
[9]^b) shows a similar comparative study at a larger value of Ri n j, closer to Ri n j,2- In this 
figure we see a distinct difference between the spectra. The additional third order term 
is apparently responsible for the splitting of the A, P, and S families, illustrated earlier 
inFig.gb). 

In Fig. [lOj we plot k\ against Ri n j(— Ri n j,i)- A linear dependence is evident. The 
slope of the line is approximately —1/3. The flow is unconditionally linearly stable above 
the line and conditionally unstable otherwise. For small values of Ri n j, we have seen in 
Fig. [T] that the principal effect is to skew the velocity profile towards the upper wall. A 
similar asymmetric skewing of the velocity profile is also induced in an annular Couette- 
Poiseuille (ACP) flow, through geometric means by varying the radius ratio, 97, (defined 
as the radius of the outer stationary cylinder to the radius of inner moving cylinder). ACP 
flow has been studied extensively by Sadeghi k, Higgi ns (1991)| and we superimpose their 



results on ours, in Fig. [101 The comparison is striking. We believe there are 2 features 
of Fig. 10 that are unusual and worthy of note. Unsurprising is of course the identical 
limits Ri n j = = (77 — 1). Note that R{ n j — > is the PCP flow, and 77 — > 1 represents 
the narrow gap limit of ACP, which is also the PCP flow. 

The first feature is the very similar linear decay in critical k = ki(Ri n j), from the PCP 
values. It can be argued along the lines of Mott fc Joseph (1968)| that for a fixed Couette 
component (k), increasing the cross-flow for the PCP flow, or the radius ratio in the ACP 
flow of Sadeghi & Higgins (1991) skews the velocity profile more towards the moving 
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Figure 9: Eigenspectrum for (k, a, R) — (0.5, 0.2, 31656) (a) Rinj — 0.5 (Critical Condi- 
tions) and (b) Ri n j — 23.5. Symbol o indicates the eigenspectrum from the O-S equation 
while □ indicates the spectrum obtained by neglecting the additional cross-flow inertial 
term. 




0.2 



0.4 0.6 0.8 1 1.2 

Rinj ,1 5 T) 1 



1.4 



Figure 10: k\ as a function of Ri n j,i (shown by □) as well as the radius ratio, r\ (shown 
by •) in ACP flow (|Sadeghi fe Higgins (1991)]> 



boundary, thus increasing asymmetry and thereby stability. Since it has been already 
observed in Fig. [9ja) that for small Rinj the influence of injection on the eigenspectrum 
is through the velocity profile only, we do expect stabilisation. However, when (77 — 1) 
and Ri n j are of 0(1), we can see no obvious quantitative relation between these flows 
and even the stability operators are quite different. 



The second noteworthy feature of Fig. 10 is that there is a minimum value of k\ (ki jm in) 



below which it is not possible to produce unconditional stability by applying (modest) 
cross-flow. This minimum value is found when k\ — > c rjCr i t . We have found approximately 
that k\^ m in — 0.19 and the corresponding Rm 1 — 1.29. This is very similar to ISadeghi 
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|fc Higgins (1991)] who found that the critical layer near the moving wall of ACP flows 
remained up to c r)Cr jt ~ 0.18. 

3.2.1. Linear energy budget considerations 



The strong analogy with the ACP results of Sadeghi & Higgins (1991) suggests that 



a similar mechanism may be responsible for the stabilisation and cut-off behaviour. To 
investigate this we examine the linear energy equation, derived in modal form from the 
Reynolds-Orr energy equation. This yields the following two identities: 



(((j) r D(j)i - <f>iD<p r )Du) 



1 

aR 



[Li + 2a 2 If + a 4 !'} 



Lf + a 2 L 2 

<(a 2 |0| 2 + \D(j)\ 2 )u) + ^<a 2 (0 r L>^ - fcDfo) + (Dcf> r D< 
_ art 

L 2 + a 2 L 2 

where L^ — ^fc(0) is the semi-norm defined by: 



(3.2) 



D<t>iD 2 ct> r )) 



(3.3) 



h = 



dy 



1/2 



and where 



f(y) dy. 



Before proceeding further, we note that Ri n j only appears indirectly in (3.2), reinforcing 
the assertion that for order unity R{ n j, the principle contribution to stability of injection 
is via the mean flow. Indeed, in the long wavelength limits of cut-off k that we have 
studied, we have found values A = (aR) -1 < 10~ 4 for instability. Thus, in (3.3) the term 
directly involving Ri n j has minimal effect on c r , explaining the observations in Fig. [9^a). 
The identity (3.2) can also be interpreted as an energy equation, in form: 



d 
dt 



(Ti) = <r 2 ) - - <t 3 ) 



R 



where 



7i 

T 3 



0.5 {\D(/)l 
O.barDu, 
0.5{\D 2 (j)\ 



dt 



Ti = aCiTi, 



2 



+ 2a z \D</)\ +a 



iD4>r 
4 I / |2x 



(3.4) 

(3.5) 
(3.6) 
(3.7) 



The left-hand side of (3.4) represents the temporal variation of the spatially averaged 



(one wavelength) kinetic energy. The first term on the right-hand side of (3.4) is the 
exchange of energy between the base flow and the disturbance. The last term, (-^ (T3)), 
represents the rate of viscous dissipation. At criticality, the two terms on the right-hand 
side balances each other, but the spatial distributions of T2 and T3/R indicate where the 
energy is generated and dissipated in the channel. 

Sadeghi h Higgins (1991)| extensively utilised this linear energy approach in studying 
the effect of k on stability of ACP flow. They found that increase in the value of k— c rjCr i t 
decreases the Reynolds stress (r) near the moving wall until it becomes negative, hence 
stabilising. The critical layer near the moving wall vanishes for k > c r , cr it and as k 
increases the Reynolds stress becomes progressively negative within the critical layer at 
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Figure 11: Distribution of energy production (T2) and dissipation (7^3) terms across the 
domain corresponding to criticality at Ri n j— (a) 0, (b) 0.2, (c) 0.4 and (d) 0.6. In all the 
cases, k — 0.5. Dash-dot-dot line with symbol □ represents T2, dashed line with filled A 
represents and solid vertical line represents the location of the critical layer. 



the fixed wall, but this behavior is destabilising since the velocity gradient is negative 
there for ACP flow. 



Figures 11 ^a)-(d) examine the distribution of T 2 and T$/R for the least stable eigen- 
mode for the parameters listed in Table [TJ i.e. we fix k — 0.5 and increase Ri 

Rinj R 



up to 

inj,i ~ 0-6- The critical layer is marked with a vertical line. We observe that 
both the rate of energy transfer and the rate of viscous dissipation decrease with the 
cross-flow. Without cross-flow, T2 is positive and negative respectively in the lower (in- 
jection) and upper (suction) halves of the domain. Increasing the cross-flow decreases 
both the positive (near injection wall) and negative (near suction wall) peaks. The loca- 
tion of the critical layer also moves away from the injection wall due to the skewing of 
the velocity profile. When Ri n j ~ Ri n j,i, (T2) and (T3) not only equalize but (since (p 
has been normalised), will have magnitudes 0(a~ 1 ) since aR — constant at cut-off; (see 
also Sadeghi & Higgins (1991)) . This reduced energy budget as Ri n j ~ Rinj,i- This is 
the primary reason for the cut-off. 
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Figure 12: Neutral Stability Curves (NSCs) for (a) k = and (b) 0.18 at different Ri n j. 
The symbols indicate different values of Ri n j and are as follows: x — > R in j — 0, o — y 
Rinj — 6 in (a) and 4 in (b), □ — > Ri n j — 12 in (a) and 8 in (b) 



3.3. Summary 

For the range of small to order unity Rinj with k ^ c rjCr it, the flow instability is dom- 
inated by long wavelength perturbations. This instability mchanism exhibits a cut-off 
phenomenon characterised by a near linear boundary in the (Ri n j, fc)-plane. The initial 
cut-off mechanism is very similar to that for ACP, as studied by Sadeghi & Higgins 
(1991) combining skewing of the velocity profile, shifting of the critical layer and decay 
of the net perturbation energy. 



4. Intermediate Ri n j and short wavelength instabilities 

We now consider the range ^ k ^ c rjCr it, in which the critical layer at the upper 
wall is still present. We investigate its stability characteristics by adding cross-flow of 
intermediate strength (0 ^ Rinj % 21), avoiding for the moment the second transition. 
It is intuitive that the presence of the critical layer will affect the stability behavior. To 
verify this we have studied the two extremities of the range of k considered, i.e. k — 
(PP flow) and k — 0.18. The respective NSCs are shown in Fig. 12 It is evident that the 



presence of the critical layers render shorter wavelength modes unstable. Yet, it is also 
observed that with Ri n j in this intermediate range, the stability increases dramatically. 

We have been unable to make any advance analytically in this range of Rinj-, and 
therefore have proceeded numerically. First we note that when we have considered k > 
0.19 for the range of 1.3 < Ri n j < 21, we have found that the least stable modes are long 
wavelength modes and that these are linearly stable. Thus, k > 0.19 appears to represent 
an absolute cut-off in this range of Ri n j- 

For smaller k we have seen that the NSCs occur with wavenumbers that are 0(1) and 
apparently increasing with Ri n j- Unlike the long wavelength problem, the asymptotic 
behaviour along the branches of the NSCs is not easily treated. At fixed large R, we are 
able to compute numerically a cut-off value of k for increasing Ri n j, i.e. k — ki(Ri n j,R). 
These cut-off curves do lie below k ~ 0.19, but are not wholly independent of R, at least 
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Table 3: Cut-off values evaluated for shorter wavelength instabilities for R — 10 6 . 



o 




Figure 13: Shorter wavelength cut-off showing k\ as a function of Ri n j,i- The flow is 
linearly stable for R ^ 10 6 above the curve. The values in Table pi are marked by □. 



within the range of R up to which our numerical code is reliable, i.e. it is quite possible 
that these asymptote to a cut-off curve as R — > oo, but we cannot reliably evaluate 
this limit numerically. As an example of this numerical cut-off, (at R — 10 6 ), we have 
computed the cut-off values Ri n j t i, as listed in Tableland shown in Fig. 13 ^a). For the 
range 1.4 < R in j < 11.8, the cut-off is close to k ~ 0.19. 



Although we see that the unstable wavenumbers increase with Ri n j in Fig. 12 note 
that asymptotically as a — > oo the short wavelengths are stable. To see this, from (3.2) 
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we bound 

{(<t> r D<t>i - <f>iD<f> r )Du) < \Du\ max I I 1 ^ 0.b\Du\ max [alQ +lf/a], 
so that Ci < provided that: 

#<^^, (4.1) 

(and better bounds are certainly possible). In Table [3] we note that the maximal critical 
wavenumber is in fact attained at an intermediate Ri n j. 

4.1. Behaviour of preferred modes for intermediate R 



inj ■ 



In our preliminary results, (^2.3), we saw that at fixed values of (R, k, a), increasing the 
Rinj led to regimes of stabilisation, then destabilisation, and then finally stabilisation. 
For k ^ c rjCr it, only long wavelengths appear unstable and how the cut-off values of k and 



Rinj vary in this regime are illustrated in Fig. 10 For the lower range of k, our results are 
primarily numerical, indicating a cut-off value k ~ 0.19 for 1.3 < Rinj ^5 H-8 an d then 
with decaying cut-off k for 11.8 < Ri n j < 20.8, as illustrated in Fig. [13] Therefore, we 
have linear stability as we cross some cut-off frontier, k > k\{Ri n j) in the (Ri n j, fc)-plane, 
(alternatively for Ri n j > Ri n j,i)- 

We now consider what happens to the certain eigenmodes (preferred modes) as we 
extend the injection cross-flow up until the second critical Rinj- Our analysis up to now 
suggests that the behaviour may be different depending on whether we consider small or 
moderate k. In Fig. fT4| we have plotted the locations of certain eigenmodes as Ri n j is 



increased, by keeping the Reynolds number R constant at 10 6 . This gives us some idea 
of how cut-off behaviour changes with Ri n j- Although the "preferred modes" are simply 
those we have selected, we implicitly mean modes that are involved in the transition from 
stable to unstable as one of our dimensionless parameters is varied (here Ri n j), i-e. at 
some point a preferred mode becomes the least stable mode and then unstable. 



Figure 14 ^a) shows two eigenmodes corresponding to k = 0, (PP flow). A least stable 
long wavelength mode is tracked for a — 0.001, denoted by 'A'. This mode is stable at 
Rinj — and its stability increases further as Ri n j increases up to around 1.7. However, 
further increases in Ri n j destabilise this mode progressively until it becomes unstable at 



Rinj — 25. In the inset of Fig. 14 ^a) we have also plotted the least stable short wavelength 
mode at a — 3.5227. Such modes become unstable only under the influence of cross-flow of 
intermediate strength. This particular mode, (denoted by 'B'), starts becoming unstable 
approximately when Ri n j > 15, but recovers stability later for Ri n i ^ 20.8. This behavior 



is a direct consequence of the trajectory of the NSCs observed in Fig. 12 ^a). The preferred 
mode 'B' is the critical mode at cut-off, (see Table p|. Thus PP flow with cross-flow is 
unconditionally linearly stable in the range 20.8 ^ Rinj ^5 25. 

For larger k, the stability behavior is primarily governed by the long wavelength modes, 
as shown in Fig. [l4|b) for k — 0.5. The least stable mode corresponding to a — 0.01 is 



unstable for Ri n j — 0, denoted mode 'C This viscous mode becomes stable when Rinj 
increases to 0.6, which is indeed the cut-off value, i.e. Ri n j,i - This is expected, according to 
Table[2] Mode 'C is weakly damped and its stability increases for Ri n j ^ 3, after which it 
starts destabilising. The mechanism of this destabilisation can probably be analysed along 



the lines of resonant interactions of the Tollmien-Schlichting (T-S) waves; see |Baines et 
al. (1996)| To show this interaction, we have traced the locus, (for Ri n j — [7,30]), of the 



least stable inviscid short wavelength mode 'D', at a — 2.5. This mode, being inviscid, 
remains stable but has q very close to zero as Ri n j increases. The wave speed c r decreases 
continuously with Ri n j for mode 'D'. The resonant interaction takes place when its wave 
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Figure 14: Behavior of preferred modes (belonging to different wavelengths and denoted 
by alphabets 'A'-'D') under the influence of cross-flow with R — 10 6 . Symbols □ and o 
respectively imply the starting and the ending position of the preferred mode in the q, c r 
plane, whereas the dots ('.') trace the locus. The difference in Ri n j between consecutive 
dots is 0.1. (a) k = 0. Mode 'A' has a = 0.001 and is traced for #^=[0,25]. Mode 'B' has 
a — 3.5227 and is traced for R in j — [15,21] (shown in the inset), the position at Ri n j — 15 
is marked by (b) k — 0.5. Mode 'C has a — 0.01 and is traced for Ri n j — [0,30]. Mode 
'D' has a — 2.5 and is traced for i^ nj -=[7,30]. 



speed matches with that of mode 'C, which signals the destabilisation of mode 'C. This 
destabilisation continues until mode 'C becomes unstable when R in j > 30. 

In Fig. [15] we show examples of the stream function for the preferred modes, corre- 
sponding to various k and Ri n j in the transitions of Fig. 14 For the long wavelength 
mode 'C, Figs. 15'a)-(c) show that strong Ri n j appears to skew the streamlines towards 
the lower wall. The same is true for the long wavelength mode 'A' under strong injection; 
T). On the contrary, Figs. 



see Fig. 



15 



15 



d)-(e) show that for large Ri n j, the streamlines 



of the shorter wavelength mode 'B' are skewed and localised towards the upper wall. 



5. Stability and instability at large Ri n j. 

We turn now to the transition to instability at Ri n j,2 and then later to stabilising effects 
at very large Ri n j. As observed in { 2.3 the transition at Ri n j,2 appears to be independent 
of streamwise Reynolds number R (see Fig. [5tb)) and occurs for all k. Although there is 
sensitivity to k, it is not very significant. Values of Ri n j,2 are found for all k G [0, 1] and 
are in a fairly tight range of R{ n j ~ 22 — 25. 

As suggested in the previous subsection, although instability at moderate Rinj may 
be either short wavelength or long wavelength, according to {k — c rjCr ^), as we approach 
Rinj,2 from below it is the long wavelengths that are unstable. Figure 16 1 shows the 
neutral stability curves corresponding to PP flow for R in j just above Ri n j,2- The NSC's 
are nested with decreasing Ri n j and as we approach Ri n j,2 the upper and lower branches 
of the NSC are seen to coalesce. The slope of the two branches suggests that a ~ R -1 
in the limit of cut-off, and hence the previous long wavelength approximation, leading 
equation (3.1), should be effective for predicting the cut-off in the (k, A)-plane; (recall 
A = (aR)^ 

Figure [l6p shows the NSC's obtained from long wave approximation. The cut-off 
velocity &2 is the maximum value of k encountered along the NSC for a given Ri n j — 
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(d) (e) (f) 

Figure 15: Isovalues of the normalised perturbation stream functions (ip) for the preferred 
modes at R = 10 6 under different Ri n j- The streamwise extent of the domain is one 
wavelength. Corresponding to k — 0.5, the long wavelength mode 'C is shown for (a) 
Rinj = 0.1 (unstable), (b) Ri n j — 1 (stable) and (c) Ri n j — 30 (unstable). Corresponding 
to k = 0, ^ for two different preferred modes, viz. 'A' and 'B' are shown. The shorter 
wavelength mode 'B' (a — 3.5227) is shown for (d) R{ n j — 15 (unstable) and (e) Ri n j — 
21 (stable). The longer wavelength mode 'A' (a — 0.001) is shown for (f) R{ n j — 25 
(unstable). 



Rinj,2- Unlike Fig. [8j the entire range of k becomes unconditionally stable. For Ri n j : \ < 
Rj n j < 22.2, Cj < V k € [0 , 1]. Another significant difference with Fig. [8] and the results 
of Cowley fc Smith ( 1985) | is that "bifurcation from infinity" is not observed as k — > 0. 



This is possibly because the curves bifurcate from infinity for negative values of k, but 
we have not studied this range. Finally, we mention that for cross-flow rates slightly 
greater than Ri n j,2, the R cr it is relatively low for the entire range of k. For example, 
Rinj,2 ~ 23.8 for k — 0.5, (implying R cr n — > oo as Ri n j — > Rj~ n j 2)- Increasing Ri n j to 25 
decreases R cr it to around 6000. Thus, on crossing Ri n j 7 2 we find a dramatic decrease in 
the flow stability. 

5.1. Linear energy balance at -Rinj, 2 
An interesting feature of transition at Ri n j,2 is the independence with respect to R. With 



reference to the energy equation (3.4), this insensitivity implies that in this range IT2I 
is much larger than the viscous dissipation, ^^3. In other words, at criticality q = 
is achieved by a balance of energy production and dissipation within T2, more so than 



via balance with the viscous dissipation. Figure 17 ^a) investigates the energy budget 



at criticality for k — 0.5 at Ri n j — 25. The critical parameters are observed to be 
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Figure 16: (a) NSC of PP flow (k = 0) when R, 
of Ri 



r: 



The different values 



vmj are 22.5 (dashed line with x), 23 (dash-dot line with □) and 24 (dash-dot- 
dot line with o). Near cut-off, aR is constant along the upper and lower branches, 
(b) Long wave NSCs showing the dependence of log 10 A on i The different values of 
Ri n j are 22.4 (□), 23.5 (o), 24 (o) and 25 (x). Cut-off is achieved over the entire range 
of k, i.e. [0,1]. 
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Figure 17: (a) Distribution of energy production (T2) and dissipation (^^3) terms across 
the domain corresponding to criticality at Ri n j — 25. Dash-dot-dot line with symbol □ 
represents T2, dashed line with filled A represents and solid vertical line represents 
the location of the critical layer, (b) Reynolds Stress r distribution at criticality for 
Rinj — (denoted by □ symbol), Ri n j — 0.6 (denoted by x) and Ri n j — 25 (denoted 
by o). The location of the critical layers are shown by solid lines with corrosponding 
symbols. 
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(a cr it, Rcrit) — (0.31, 6000). This implies that crossing the cut-off Ri n j,2, there is a tran- 
sition from unconditional stability (R cr it — > oo) to high instability (R cr it — 6000). Com- 
paring with Fig. [IT] (which shows energy distribution corresponding to criticality for 
k — 0.5 and Ri n j ^ Rinj,i), it is obvious that T2 has a higher amplitude while the 
viscous dissipation j^T 3 is weaker. 

This behavior is due to the generation of larger Reynolds stresses r, as Rinj increases, 
as illustrated in Fig. [T7|fo). The dominance of T2 over the viscous dissipation suggests 
that the critical layers have little to do with instability in this range. Note that r is small 
in the critical layer, which has now moved towards the channel centre, and hence T2 is 
also small. Referring to Fig. (2^b), the vanishing vorticity gradient (D 2 u) found in the 
bulk of the flow domain at high values of Ri n j removes/diminishes the singular effects 
associated with the critical layer. 

The growth of r is probably not responsible for the spreading of the spectrum along 
the real axis, that we have observed in Fig. [4Tb). Equation (3.3) may be rewritten as: 



{{a 2 \ct>\ 2 + \D<j>\ 2 )u) + ^(a 2 r - ^D 3 ^ + faD 3 ^) 
/? + a 2 ! 2 



(5.1) 



The first term leads simply to values of c r in the range of u. The second term does contain 
a 2 r, i.e. longitudinal gradients of the Reynolds stresses. However, note that even for the 
shorter wavelengths we have a ~ 0(1), and if we consider long wavelengths, we have 
typically found instability only for aR 3> 1. Thus, even for these larger Ri n j, the term 
involving a 2 r is likely to be insignificant. 

The extension of c r beyond the usual bounds of the base flow velocity is therefore 



due to the 3rd derivative terms in (5.1), which cannot be bounded by the denominator. 
Interestingly therefore, the larger values of c r , which indicate less regular eigenmodes, 
also lead to larger viscous dissipation, and hence more stable modes. This explains the 
shape of the spectrum in Fig. [4|b). 



In Fig. 14 ^b), we tracked the behaviour of mode 'C as Ri n j increased. This mode 
becomes unstable for Ri n j ^ Rinj, 2, implying that it governs the transition behavior. 
In Fig. 18 we show the evolution of the energy balance terms for this mode as Ri n j 



increases from zero. This mode is stable from Ri n j — 0.6 to 30. The cut-off achieved 
at Ri n j = 0.6 is primarily due to the increased viscous dissipation at both walls. This 
phenomenon continues until Ri n j ~ 3, see Fig. 18 'd). At this point, the 'viscous hump' 
observed near the lower wall gets amplified. This mechanism is probably due to the 
resonant interaction between mode 'C and an (approximately) neutrally stable inviscid 
mode, for example mode 'D'. Further increase in Ri n j thins out the viscous layer at the 
suction wall faster than that at the injection wall. Suction negates both the exchange as 
well as the dissipation of energy, and the viscous hump is localised within the lower half of 
the channel, i.e. injection side. j^T 3 reduces faster than T2 and finally the mode becomes 
unstable when Ri n j increases to 30. The condition at this point is (T2) > \ (T3); see 



Fig. 18 [i). It is interesting to observe that the mode becomes unstable when the viscous 
hump reaches the centre of the channel. Further increase of Ri n j results in a gradual 
reduction of T2 and the mode becomes stable again. The mean perturbation kinetic 
energy q(y) distribution provides further insight into the instability mechanism. It is 
defined modally to be: 

g=i(W| 2 + a 2 H 2 ) (5.2) 
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Figure 18: Distribution of energy production (T2) and dissipation (7^3) terms across the 
domain corresponding to mode 'C at Ri n j — (a) 0, (b) 0.6, (c) 1, (d) 3, (e) 10 and (f) 30. 
Dash-dot-dot line represents T2, solid line represents j=>?3. 



Figure 19 shows the mean perturbation kinetic energy profiles for mode 'C at different 
Rinj- Each distribution of q has been normalised by its maximum value. 

Without any cross-flow, the amount of energy in the two halves of the domain are 
comparable, the suction half having ~ 43% of the energy, (note k — 0.5). Increasing 
cross-flow up to Ri n j ~ Rinj, 1 = 0.6 increases the secondary peak until the cut-off is 
achieved. The energy in the suction half at this point is 46.7%. The primary peak moves 
toward the lower wall but cannot reach it because of the no-slip conditions. At Rinj > 1, 
the primary peak starts moving away from the suction wall. At Rinj — 3, the mode 
is at its maximum stability (see Fig. [l4^b)) . At this point, the perturbation energy is 
highly localised within the lower |th of the channel, along with a small secondary peak 
at the upper quarter. Further increase of Rinj to 10 causes the secondary peak to vanish; 
the energy content in the suction half being only ~ 7.6%. The resonant interactions of 
T-S waves result in the development of a secondary peak from the primary peak itself. 
During this process, the secondary peak slowly separates from the primary peak and 
moves in the direction of the upper wall. For R in j — 30, the perturbation reaches the 
channel centre and the mode becomes unstable. The amount of energy in the suction 
half increases to 18.1%. For even higher values of Rinj, for example 45, the upper half 
holds ~ 32% of the energy. 

Thus it appears that the onset of the cut-off at Ri n j,i occurs when the secondary peak 
holds maximum energy. Increasing injection decays this peak until it reaches a minimum 
and then starts to grow out from the primary peak. The end of the cut-off regime, marked 
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Figure 19: Non-dimensional mean perturbation kinetic energy profiles for mode 'C at 
different Ri n j- Solid lines with symbols denote the unstable modes. For each Ri n j, q has 
been scaled by its maximum value. 



by Ri n j > Rinj, 2, occurs when the secondary peak reaches the channel centre and holds 
sufficient energy. 



5.2. Eventual stabilisation at Ri n j,3 

We have not studied in detail the eventual stabilisation of the flow at very large Ri n j (i-e. 
Rinj ~ -Rmj, 3) 5 but we believe the energetics of this stabilisation are due to a decay in the 



energy production. This can be seen most clearly from the identity (3.2), which is in the 
same form as that for any parallel shear flow, i.e. cross-flow only influences (3.3) directly. 
Joseph has used this expression to derive general bounds that depend on \Du\ max , and 
various functional inequalities; see Joseph (1968), Joseph (1969) For example, we have 
linear stability provided that: 



aR\Du\max < max(£ 2 7r + 2 3/2 a 3 , £ 2 7r + a 2 7r) 



(5.3) 



where £ = 2.36502 is the least eigenvalue of a vibrating rod with clamped ends at y = ±1. 



The condition (5.3) evidently holds for the flows we consider, but is very conservative 

This conservatism at large Ri n j stems directly 



and especially so in the limit of large Ri n j 
from the simplistic treatment of Du in bounding the energy production term 

(((f) r D(f)i - (j)iD(f) r )Du) < \Du\ m axhh- 



With reference to Figs. [I] Sz [2] and to (2.11), we see that at large Ri n j the base velocity 
profile consists of a thin layer near the upper suction wall, within which Du ~ \Du\ ma x ~ 
Rinj, which has thickness of 0{R~ n 1 -). Away from this thin boundary layer, the velocity 
gradients are of size Du ~ 2(kRi n j)~ 1 + 0(Ri n je~ Rinj ( 1_y )) . Note however, that within 
this suction layer, we have ~ (1 — y) 2 due to the boundary conditions on the pertur- 
bation. Therefore, taking a nominal suction layer boundary at y — y s , we may estimate 
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Figure 20: Variation of k with Ri n j. The filled □ symbols show the long wavelength cut- 
off achieved for 0.7 ^ k ^ 0.19. The filled o symbols show the shorter wavelength cut-off 
for 0.19 > k ^ evaluated numerically for R — 10 6 . The filled o symbols imply the second 
long wavelength cut-off. The shaded region depicts the entire zone of unconditional linear 
stability. 



as follows: 

(((f) r D(f)i - (j)iD(j) r )Du) = 



i-y s 

l 
2 



kR 



r D(j)i 
•l 



<t>iD(t> r )Du dy + 
,D(f>i 



l-Vs 

Dcj) r \ dy + 0(\Du\ max 



kR 



Ioh + 0{Rr*) 



i nj 



Following Joseph (1969)| this leads directly to the bound 

2aR 



kR, 



< max(£ 2 7r + 2 3/2 a 3 , ^ 2 tt + a 2 7r) 



r Dcf)i — (j>iD(/) r )Du dy 

(i - y s f) 

(5.4) 
(5.5) 



sufficient for linear stability at large Ri n j, (with asymptotically kR{ n j > 4 required). In 
other words, at large Ri n j, the energy production T2 will decay like (kRinj) -1 at leading 
order, so that the viscous dissipation need only be of this order to stabilise the flow. 



6. Summary 

To summarise, we have presented a detailed analysis of linear stability and instability in 
the (Ri n j , k) -plane, for PCP flow with cross-flow. The most complete analysis concerns 
the important range of low Ri n j and modest k. In this range we have demonstrated 
that the stabilisation mechanism, due to either injection or wall motion, is essentially 
the same. Long wavelengths dominate. Skewing of the velocity profile shifts the critical 
layer and at the same time the energy production is diminished until viscous dissipation 
dominates at cut-off. In Fig. flO] we have also shown an interesting quantitative analogy 



with the cut-off behavior of ACP flows; see Sadeghi & Higgins (1991) 

This lower range of Rinj and modest k is probably that which is most important 
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practically. Essentially, this range allows one to compensate cross-flow by wall-motion and 
vice-versa, achieving unconditional linear stability via either mechanism. With reference 
to Fig.[TJ it is the range of Ri n j in which the cross- flow and wall motion are modifications 
of a base Poiseuille flow. Due to the scaling, the peak velocity is always 1, but at larger 
Ri n j with modest k the Poiseuille component is completely dominated by cross-flow and 
wall motion. 



Globally, the cut-off regimes in the (Ri n j, &)-plane are as illustrated in Fig. 20 The 
shaded area shows the region of unconditional linear stability. In the intermediate range of 
approximately 1.3 ^ Rinj ^20.8 values of k > 0.19 are dominated by long wavelengths 
and are stable. Below this value, we are able to compute numerical cut-off curves for 
fixed R. With the limits of our computations, we cannot determine if these cut-off curves 
asymptote to an unconditional cut-off curve as R — > oo. 

There appears to be a short band of unconditional linear stability for all computed 
values of k around approximately 20.8 ^ Rinj ^ 22, before the destabilisation occurs at 
larger Ri n j — Rinj, 2- Since this band can make PP flow unconditionally stable, it could be 
effectively used in applications where wall motion is not feasible, e.g. cross-flow filtration, 
medical dialysis. From the practical perspective, it is worth noting that the transition 
across Ri n j,2, is from unconditional stability to critical values of R which are relatively 
modest (e.g. in the range 10 3 — 10 4 ) just a short distance beyond Ri n j,2- Assuming that 
the PP flow is linearly unstable, this means stabilisation can be achieved with cross-flow 
velocities of the order of 1% of the mean axial flow velocity. 

This destabilisation at Ri n j,2 is again a long wavelength mechanism, which we have 
analysed using the long wavelength approximation of Cowley Sz Smith (1985)| A possible 



cause of this instability has been found to be resonant interactions of the T-S waves. 
Study of the linear energetics of the upper limit, Ri n j, 2, has shown that neither viscous 
dissipation, nor the involvement of a critical layer are significant. Rather, the balance 
of energy production and dissipation within T2 keeps the mode neutrally stable. Energy 
analysis of the preferred mode 'C has revealed that the precursor of the transition 
to instability from unconditional stability is the amplification of disturbances near the 
injection wall. The mean perturbation kinetic energy has also been analyzed. It has shown 
that the lower limit occurs when the secondary peak holds maximum energy. Increasing 
injection decreases the secondary peak until it reaches a minimum and then it starts to 
grow from the primary peak. When the secondary peak reaches the channel centre and 
holds a sufficient amount of energy, the unconditional stability mechanism breaks down. 

The final stabilisation occurring at large Ri n j ^ Rinj,3 has been analyzed using linear 
energy bounds. By careful treatment of the energy production term, we are able to show 
that the energy production terms decreases asymptotically like RJ^j as Rinj — > 00. We 
believe that this mechanism leads to the eventual domination of the viscous dissipation 
at large enough Ri n j- 

In terms of the spatial structure of the perturbations, we note that the stabilisation 
at small and moderate Ri n j are both long wavelength phenomena for which the approx- 
imation of Cowley Smith has been shown effective. Implicitly therefore, the critical 
wavenumbers scale like R~ x in these limits. For the shorter wavelength instabilities we 
have not analysed the asymptotic behaviour of the wavenumber with R. A more detailed 
look at the spatial structure of certain eigenmodes has been presented in Fig. [15} This 
shows a skewing of the streamline recirculatory regimes towards the lower wall for long 
wavelengths as Ri n j is increased, and towards the upper wall at shorter wavelengths as 
Rinj is increased. 
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